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DEFORMATION OF QUADRILATERALS AND ADDITION 
ON ELLIPTIC CURVES 

IVAN IZMESTIEV 


Abstract. The space of quadrilaterals with fixed side lengths is an 
elliptic curve. Darboux used this to prove a porism on foldings. 

In this article, the space of oriented quadrilaterals is studied on the 
base of biquadratic equations between their angles. The space of non- 
oriented quadrilaterals is also an elliptic curve, doubly covered by the 
previous one, and is described by a biquadratic relation between diag¬ 
onals. The spaces of non-oriented quadrilaterals with the side lengths 
(oi, Gt 2 , (I3, 04 ) and (s — 01, s — 02 , s — 03, s — 04 ) turn out to be isomorphic 
via identification of two quadrilaterals with the same diagonal lengths. 

We prove a periodicity condition for foldings, similar to Cayley’s 
condition for the Poncelet porism. 

Some applications to kinematics and geometry are presented. 


1. Introduction 

1.1. The Darboux porism. Let ABCD be a planar quadrilateral. De¬ 
note by D' the image of D under reflection in the line AC. We call the 
transformation 

F d : ABCD ha ABCD' 

the D-folding. It acts on the set of all quadrilaterals ABCD with A / C. 
Define the C-folding in a similar way, see Figure [l] 



Figure 1 . Composition of C- and D-foldings. 

Alternate the C- and the D-foldings, so that the vertices A and B stay 
in their places, while C and D jump. If we are lucky, then after several 
iterations the points C and D return to their initial positions. 
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Definition 1.1. A quadrilateral ABCD is called n-periodic, if it is invari¬ 
ant under the n-fold iteration of the composition of C- and D-foldings: 

(. F c o F D ) n (ABCD) = ABCD 

The quadrilateral from Figure [l] is 3-periodic, see Figure [2] Surprizingly 
enough, the periodicity depends only on the side lengths. 




Figure 2. A 3-periodic quadrilateral. Side lengths are l,3,3\/5, 5. 

Theorem 1 (Darboux |9j|). If a quadrilateral is n-periodic, then every 
quadrilateral with the same side lengths is n-periodic. 

The clue to the Darboux porism is the fact that the space of congruence 
classes of quadrilaterals with fixed side lengths is an elliptic curve, and the 
foldings Fq and Fjj act on it as involutions. Hence Fc o Fjj is a translation, 
just as well as (FcoFc) n - If a translation has a fixed point (one quadrilateral 
is n-periodic), then it is an identity (all quadrilaterals with the same side 
lengths are periodic). 

Unlike the Poncelet porism im, its relative, the Darboux porism is much 
less known and was rediscovered several times in the recent decades. 

In the present article we take a closer look at the configuration space 
Qo r (a) of quadrilaterals with side lengths bGK 4 . We study also the space 
Q c (a ) of non-oriented quadrilaterals, which is an elliptic curve doubly cov¬ 
ered by Q^ r (a). There are some surprising algebraic identities that result in 
a natural isomorphism between the spaces Q c (a) and Q c (a), where a^ais 
a certain involution on M 4 . Geometrically this leads to “conjugate pairs” of 
quadrilaterals. Finally, we express n-periodicity of a quadrilateral in terms 
of its side lengths. 

1.2. Euler-Chasles correspondence and Jacobi elliptic functions. 

Introduce the variables 

. <Pi . ¥>2 

z\ = tan —, 22 = tan — 
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where p\ and <^> 2 are adjacent angles of a quadrilateral. Then the space of 
congruence classes of quadrilaterals, respecting the orientation, with fixed 
side lengths becomes identified with the algebraic curve 

(1) C22^1~2 + c 20 z l + c 02 z 2 + ^ c ll z l z 2 + C()0 = 0 

Here coefficients Cij depend on the side lengths. This is a special biquadratic 
equation in two variables. Biquadratic equations are also known under the 
name of Euler-Chasles correspondences [7j. In the past decades they at¬ 
tracted a lot of attention: as a basis for QRT- maps in the theory of discrete 
integrable systems [12]; as a solution of the Yang-Baxter equation [UG±2|; as 
an approach to flexible polyhedra nana ng . 

It is known that in a generic case the curve ([!]) can be parametrized as 

Zl=pi<p(t), Z 2 =P 2 V(t + T) 

where ip is a second-order elliptic function with simple poles and zeros. Since 
any pair of adjacent angles is related by a biquadratic equation, all angles 
turn out to be scaled shifts of the same elliptic function. 

If a = (oi, 02 , 03 , a 4 ) are the side lengths of a quadrilateral, then the 
non-degeneracy condition ensuring that the configuration space Q 0 r(a) is 
an elliptic curve is 

CL\ ± 02 A 0,3 i 04 7 ^ 0 

for all choices of signs. Put differently, this means o m i n + a max / s , where 
a m i n and o max are the lengths of the shortest and the longest side, and 
s = ai+ffl 2 jj a3+ffl4 is the half-perimeter. 

In kinematics, the inequality a m i n + a max < s is knows as the Grashof 
condition. It means that the shortest side can make a full turn with respect 
to each of its neighbors; also it means that the real part of the config¬ 
uration space Qor ( 0 ) has two components, that is the quadrilateral can¬ 
not be deformed into its mirror image. The following theorem describes 
a parametrization of Q^ r (a) and underlines the difference between Grashof 
and non-Grashof quadrilaterals. 

Theorem 2. If a± ± a 2 ± <23 ± gq 7 ^ 0, then the complexified configuration 
space Q^ r (a) of quadrilaterals with side lengths a is an elliptic curve C/A. 

(1) If a m i n + a max < s, then the lattice A is rectangular, and the cotan¬ 
gents of the halves of the exterior angles can be parametrized as 

cot ^ = qi dn(t + ti;k), cot ^ = <72 dn(f + f 2 ; k), 

cot = <73 dn(t + t;s] k), cot y = <24 dn(t + tp, k) 

Here k = amplitudes qi are real or purely imaginary, 

and the shifts t t satisfy 

Im ti G {0,2 1\'}, t\ - f 3 ,t 2 ~ t 4 e {±K} 
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The real part of the configuration space consists of two components 
corresponding to Im t E {0,2 K 1 }. 

(2) If Q min + fl mav > s, then the lattice A is rhombic, and the cotangents 
of the halves of the exterior angles can be parametrized as 

cot -y- = qi cn (t + ti; k), cot ^ = g 2 cn(f + t 2 ; &), 

cot ^ = g 3 cn(t + f 3 ; fc), cot ^ = g 4 cn(f + f 4 ; fc) 

LTere k = eac/i o/ i/ie amplitudes qi is real or purely 

imaginary, and the shifts ti satisfy 

Im ti E {0, 2 K'}, t\ — t 3 , f 2 — i 4 E {± K} 


The real part of the configuration space has one component corre¬ 
sponding to t E M. 


The proof of Theorem [2] is contained in Sections 3.3.4 and 3.3.5 There 
one can also find exact values of the amplitudes and shifts. The numbers a t 
appearing in the formula for the Jacobi modulus are defined as follows. 


( 2 ) 


a i = 


—ai + a 2 f ®3 f fl 4 


= s — a\, ai = s — ai , 


2,3,4 


The pair of opposite angles </?i, </? 3 in a quadrilateral is also subject to a 
relation of the form ([Tj) , but with a vanishing coefficient at z\zz- A scaling 
of variables brings the equation into the form 

(3) u 2 + v 2 = 1 + mu 2 v 2 


with a real m < 1. This curve can be parametrized as 


u = f(t-,k), v = f(t + K;k) 

where / = dn if m > 0 and / = cn if m < 0. Equation Q is closely related 
to the equation 

2 I 2 2 , 2 2 2 

u +v = a + a u v 

suggested by Edwards jl3] as a new normal form for elliptic curves. Edwards 
constructs a parametrization of u and v “from the scratch” using a version 
of theta-functions. 


1.3. The space of non-oriented quadrilaterals. Let Q c (a) be the space 
of congruence classes, disregarding the orientation, of quadrilaterals with 
side lengths a. Clearly, the map Q^ r {a) —> Q c (a) is a double cover, easy 
to describe in terms of the holomorphic parameter of Theorem [2j There is 
an unexpected natural isomorphism between the spaces Q c (a) and Q c (a) 
with a as in for every quadrilateral with the side lengths a there is a 
quadrilateral with the same diagonal lengths and the side lengths a. 
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Theorem 3. J/ai±a 2 ±a 3 ±a 4 7 ^ 0, then the complexified configuration space 
Q c (a) of non-oriented quadrilaterals with the side lengths a is an elliptic 
curve C/A with a rectangular lattice A. 

There is a natural isomorphism 

Q c (a) Q c (a) 

that identifies two quadrilaterals with the same diagonal lengths. 

If Qg r (a) has a rectangular lattice, then Qf yr (a) has a rhombic lattice. 
In particular, the above isomorphism does not lift to the spaces of oriented 
quadrilaterals. The double covers 

Qor( a ) Q C ( a ) — Q C (d) < 5 — Qor(h) 

look as shown on Figure [3| 

Theorem [3] is proved in Section |4j 



Figure 3. Spaces Q^ r (a) and Q^ r (a) covering Q c (a) = Q c (a). 


The isomorphism Q c (a ) = Q c (d) turns out to be equivalent to the Ivory 
theorem, see Figure [12} Theorem [3] holds also for spherical and hyperbolic 
quadrilaterals. This extends the Ivory theorem to the sphere and the hyper¬ 
bolic plane, with ellipses and hyperbolas defined geodesically, see Theorem 

m 

Computations that allow to express in a particularly nice way the Jacobi 
modulus, amplitudes, and shifts in Theorem[2j and those leading to Theorem 
[3] are based on a number of identities for dual quadruples (a, 6 , c, d) and 
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(a, b , c, d). These are collected below. 

ab — cd = (s — a — c)(s — b — c), ab — ab = s(s — c — d) 
abed — abed = s(s — a — 6 )(s — b — c)(s — a — c) 
ab + cd = ab + cd, a 2 + b 2 + c 2 + d 2 = a 2 + b 2 + c 2 + d 2 

ab — cd = - (c 2 + d 2 — a 2 — b 2 ) 

1.4. The periodicity condition. Benoit and Hulin |5] studied the period¬ 
icity condition by constructing a pair of circles whose Poncelet dynamics is 
equivalent to the folding dynamics of the quadrilateral. 

The following t heor em deals with the periodicity disregarding the orienta¬ 
tion. Proposition 5.4 describes how the period lengths on the curves Q c (o) 
and Qo r (a) are related. 

Theorem 4. A quadrilateral with the side lengths a is n-periodic disregard¬ 
ing the orientation if and only if the following condition is satisfied. 


A 2 

A 3 

A 3 

A 4 

Ak+i 

Ak +2 ■ ■ • 

A 3 

a 4 

A\ 

a 5 

Ak+i 

Ak +2 ■ ■ • 


Ak+i 

Ak+2 

Aik 

Ak+i 

Ak+2 

A2k -1 


= 0 , if n = 2 k + 1 


= 0 , 


if n = 2k 


Here Ai are the coefficients of the expansion 

OO 

\J {x — A (a) 2 )(x — 6 (a) 2 )(x — S(a) 2 ) = ^ AiX 1 

where 


i =o 


A (a) — a±a^ + q. 2 ® 4 ) 5(o) — a±a^ — ci 2®4 


Note that A (a) = A(a) due to the identities from Section 1.3, so that the 
quadrilaterals with the side lengths a have the same folding period (disre¬ 
garding the orientation) as the quadrilaterals with the side lengths a. This 
is also obvious from the isomorphism Q c (a) = Q c (a ) via equal diagonal 
lengths. 


1.5. Acknowledgments. Parts of this work were done during the author’s 
visits to IHP Paris and to the Penn State University. The author thanks both 
institutions for hospitality. Also he wishes to thank Arseniy Akopyan, Udo 
Hertrich-Jeromin, Boris Springborn, and Yuri Suris for useful discussions. 
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2. The space of oriented quadrilaterals in terms of their 

ANGLES 

2.1. Notation. A planar quadrilateral is for us an ordered quadruple of 
points ( A , B , C, D ) in the euclidean plane such that 

A^B, B ytz C, C^D, D^A 

Two quadrilaterals ABCD and A'B'C'D' are called directly congruent if 
there exists an orientation-preserving isometry /: M 2 —> R 2 such that 

f(A) = A r , f(B) = B', f(C) = C', f(D) = D' 

In this article, we study the set of (direct) congruence classes of quadri¬ 
laterals with fixed side lengths. A mechanical interpretation of this is the 
configuration space of a four-bar linkage with the positions of two adjacent 
joints fixed and whose bars are allowed to cross each other. 

For a quadruple a := (ai, gq, 03 , gq) of real numbers there exists a quadri¬ 
lateral with side lengths a if and only if the following inequalities hold: 

(4a) ai> 0 Vi 


(4b) ai < aj + + ai Vi, 

where {j. k, 1} = { 1 , 2 ,3,4} \ {i} in the second line. 


Definition 2.1. Fora 6 M 4 satisfying the conditions (4a) and (4b), denote 
by Qor (o) the set of direct congruence classes of quadrilaterals ABCD with 


AB = ai, BC = 02 , CD = 03 , DA = 04 


Remark 2.2. We require a congruence to preserve the marking of the ver¬ 
tices (or, equivalently, the marking of the sides). This does matter only if 
the sequence a is symmetric under the action of an element of the dihedral 
group on (1, 2, 3, 4). 

Denote by ipi the angle between the sides marked by i — 1 and i(mod 4). 
More exactly, <pi are the turning angles for the velocity vector of a point that 
runs along the perimeter in the direction given by the cyclic order (1, 2, 3,4), 
see Figure |4j left. 





Figure 4. Notations; finding a relation for opposite angles; 
finding a relation for adjacent angles. 
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2.2. Equations relating the angles of a quadrilateral. 

Lemma 2.3. The cosines of opposite angles of a quadrilateral are subject 
to a linear dependence, with coefficients depending on the side lengths. 

Proof. By expressing the diagonal length y on Figure [4] with the help of the 
cosine law first through a\, a 2 , <pi and then through < 23 , < 24 , <^ 3 , we obtain 

( 5 ) 04 + 04 + 20404 cos <p\ = a 2 + a 3 + 20203 cos +3 = 0 

□ 


This result is classical. Bricard | 6 | mentions it as well-known. The fol¬ 
lowing substitution also appears in | 6 j. 

( 6 ) Zi := tan y 

Proposition 2.4. The tangents of the opposite half-angles of a quadrilateral 
satisfy the following algebraic relation 

(7) b 22 z\z\ + b 20 zj + b 02 zl + 6 00 = 0, where 


bo 2 = (o 1 + 02 — 03 — 04) (04 — 02 + 03 — 04) 

b 2 0 = (0.1 + 02 + 03 — 04) (04 — 02 — 03 — 04) 

bo 2 = (q 1 — 02 + 03 + 04) (04 + 02 — 03 + 04) 

boo = ( a l — 02 — 03 + 04) (04 + 02 + 03 + 04) 


1 — Z z 

Proof. Substitute cos <pi = ——\ and cos +3 = 


tation yields Q. 


1+2T 


1 —z\ 
1 -\-z. 


2 

I in 


©• A simple compu- 

□ 


With a bit more work one finds a relation between pairs of adjacent angles. 


Proposition 2.5. The tangents of the adjacent half-angles of a quadrilateral 
satisfy the equation 

(8) C 22 z\zl + C2C+1 + C 02 Z 2 + 2 cnzi ^2 + coo = 0, where 

C22 = (fll — 02 — 03 — O4) (01 — 02 + 03 — 04) 

C20 = (ai + 02 + 03 — O4) (01 + 02 — 03 — 04) 

C02 = (ai — 02 + 03 + O4) (01 — 02 — 03 + 04) 

Cll = —40204 

c oo = (°i + 02 — 03 + af)(ai + 02 + 03 + 04) 

Proof. Denote by the vector running along the +th side of the quadrilat¬ 
eral, see Figured right. Expanding the scalar product on the left hand side 
of ||u4 + ui + v 2 If = ||u3|p we obtain 

(04 + 02 — a| + a|) + 20404 cos +4 + 20402 cos (p2 + 20204 cos(</?4 + ip2 ) = 0 
The substitution and a tedious computation produce ([8]) . □ 
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2.3. Bihomogeneous equations and algebraic curves in CP 1 x CP 1 . 

The substitution @ identifies M/27r, which is the range of pi, with MP 1 , 
which is the range of Z{. Thus it is geometrically reasonable to consider 
equations 0 and ([8]) as equations on (MP 1 ) 2 . The proper setting for this 
are bihomogeneous polynomials. That is, we introduce projective variables 
(zi : wi) and (z2 : W2) and rewrite equation ([8]) as 

(9) c 2 2 z\z\ + c 2 ozfw% + c 0 2 iof z\ + 2cn^iu;iZ2W , 2 + coqw\wI = 0 

This point of view doesn’t affect the affine part of the curve but may change 
the number of points at infinity (the infinity of MP 1 x MP 1 is the union 
of two projective lines instead of one for MP 2 ). Indeed, while the usual 
projectivization of ([8]) has two points (1:0:0) and (0 : 1 : 0) at infinity, 
the curve 0 has four: 

N 

For a generic choice of coefficients c^j the curve Q is non-singular, as op¬ 
posed to the usual projectivization of Q. 



2.4. A system of six equations. By Propositions 2.4 and 2.5 the angles 
of every quadrilateral satisfy a system of six equations: two of the form ([T]) 
and four of the form ([8]). In this section we show that, vice versa, under a 
certain genericity assumption every solution of the system corresponds to a 
quadrilateral. 


Proposition 2.6. Assume that the quadruple a is not made of two pairs of 
equal adjacent numbers: 

(10) neither a± = 02, <23 = 04 nor ai = aq, <22 = 03 

Then every solution (z\, Z2, Z3, Z4) of the system of six equations on the pairs 
of angles corresponds to a unique quadrilateral in Q or (a). 


Proof. By reverting the argument in the proof of Proposition 2.5 
that for every solution of 


one sees 
(a) with 


there is a unique quadrilateral in Q, 
angles ip\ and (p2- The same is true for the other three equations relating 
adjacent angles. For every solution of ([7]) there is also a quadrilateral with 
angles p\ and p 3. This quadrilateral can be non-unique only if the second 
and the fourth vertices coincide, for which a\ = 04 and 02 = <23 is needed. 

Thus for every solution (21, Z2, z 3, Z4) of the system of six equations there 
are six quadrilaterals Cij G Qor(a) with the property that in C t j the angles 
Pi and ipj have the correct values. Our goal is to show that there is C with 
all correct angles. Consider 6/2, C13, and C14. They all have the same p\. 
Assumption (10) implies that there are at most two quadrilaterals in Q or (a) 
with a given p\. Thus at least two of the three quadrilaterals must coincide. 
This yields a quadrilateral C\ with three correct angles, one of which is p\. If 
the fourth angle pi is incorrect, then again, consider the three quadrilaterals 
where this angle is correct and obtain a quadrilateral Ci with three correct 
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angles, one of which is <pi. As C± and C* have two angles in common, they 
coincide, and C = C\ = Q is the desired quadrilateral. □ 

2 . 5 . Birational equivalence between ([ 7 ]) and Q. We have just seen 
that the configuration space Q or {a) is an algebraic curve in (MP 1 ) 4 given by 
six equations. It will often be convenient to restrict our attention to a single 
equation in two variables. The following statement allows us to do so. 

Proposition 2.7. The projection of Q or (a) C (MP 1 ) 4 to every coordinate 
plane ( Zi,Zi + \) is a birational equivalence. 

The projection of Q or (a) to the coordinate plane (zi,Zi + 2) is a birational 
equivalence unless at = ai-± and aj+i = aj+ 2 . (All indices are taken in 
{ 1 , 2 , 3 , 4 } modulo A.) 

Proof. It suffices to consider the case i = 1 . As noted before, the projection 
of Qor (o) to the (zi, z 2 )-plane is injective. It suffices to show that for every 
(zi, Z2, Z3, Z4) € Qor (a) the values of z 3 and Z4 are rational functions of zi 
and Z2 ■ By projecting the quadrilateral to its second side and to the line 
orthogonal to it, we obtain 

ai cos ip 2 + a 2 + a 3 cos (p 3 + a 4 cos(y> 4 + (p 2 ) = 0 


a 1 sin tf 2 — a 3 sin (p 3 + a 4 sin(<^ 4 + (p 2 ) = 0 


It follows that 


23 = 


1 - COS g> 3 a\ COS If2 + 02 + <23 + cos((/9i + (p 2 ) 


01 sin (p 2 + a 4 sin(<^i + ip 2 ) 


sm ip 3 

which is rational in z\ and z 2 - 

Similarly, if either a\ / a 4 or a 2 / o 3 , then the projection of Q or (a ) to the 
(z\ , z 3 )-plane is injective. To show that the inverse map is rational, rewrite 
( fUj ) as 

(ai + a 4 cos (pi) cos <p 2 — a 4 sin ipi ■ sin ip 2 = — a 2 — 03 cos </? 3 
a 4 sin (pi ■ cos (p 2 + (a 4 + a 4 cos (pi) sin (p 2 = a 3 sin </? 3 

This is a system of linear equations on cos (p 2 and sin p 2 with the determinant 
(a 4 +a 4 cos (pi ) 2 +a 4 sin 2 (pi. The determinant vanishes only for pi = i: and 
a 4 = a 4 , which also implies a 2 = a 3 . If it does not vanish, then by solving 
the system we can express cos ip 2 and sin (p 2 as rational functions of zi and 
z 3 . □ 

Now that the algebraic structure of Q or (a) became clear, we can define 
the complexified configuration space. 

Definition 2 . 8 . The complexified configuration space Qg r (a) of oriented 
quadrilaterals with edge lengths a is the algebraic curve in (CP 1 ) 4 defined by 
the six equations of the form 0 and ([8]). 

Qo,. (a) is birationally equivalent to the curve ^ 


2.7 


Due to Proposition 
and is birationally equivalent to 0 unless ai = a 4 and a 2 = a 3 . 
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3 . PARAMETRIZATIONS OF THE CONFIGURATION SPACES OF ORIENTED 

QUADRILATERALS 

3 . 1 . Classifying linkages by their degree of degeneracy. The shape 
of the configuration space Q^ r (a) turns out to depend on the number of 
solutions of the equation 


( 12 ) 


a\ A <22 i i 0,4 = 0 


It is easy to see that if (12 ) has at least two solutions, then (04,02, a 3,04) 
consists of two pairs of equal numbers. This leads to the following classifi¬ 
cation of quadrilaterals. 

Definition 3 . 1 . A quadrilateral with side lengths 04, 02, (13, 04 (in this cyclic 
order) is said to be 

• of elliptic type, if equation (12) has no solution; 

• of conic type, if the equation (12) has exactly one solution; 

• an isogram, if its opposite sides are equal: 04 = 03, 02 = 04; 

• a deltoid, if two pairs of adjacent sides are equal: «i = 02,03 = <24 
or 04 = 04, <3.2 = 03; 

• a rhombus, if all sides are equal. 

Let us look at the last three simple cases before we proceed to the more 
interesting conic and elliptic types. 


3 . 1 . 1 . The rhombus. Equation 
the bihomogeneous form) 


becomes zj = z\, and Q becomes (in 


WiW 2 (ziZ2 ~ UJ1IV2) = 0 

The other equations can be obtained by cyclically permuting the indices. 
The configuration space consists of three lines 

z\ = 00, z 2 = 00, Z\Z 2 = 1 

The first two consist of “folded” configurations, when two opposite vertices 
are at the same point, and the edges rotate around this point; the third line 
corresponds to actual rhombi with ip\ + (p 2 = n. 

Note that the condition ( 10 ) is violated, and the system has “fantom” 
solutions that don’t correspond to any quadrilateral, e. g. z\ = Z 3 = 00, 
2:2 = Z4- 


3 . 1 . 2 . The deltoid. Assume a± = a 2 / 03 = 04. Equation 
(01 - 03)^1 - 2a 3 2i2 2 + (ai + a 3 ) = 0 


becomes 


Besides, in the bihomogeneous form the factor w 2 appears, so that the con¬ 
figuration space consists of two lines 

(ai - az)zl + (ai + a 3 ) 


Qm-(a) = { z 2 = 00} U < z 2 = 


2a 3 z 1 
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Again, the first line corresponds to a folded deltoid, and the second expresses 
an angle at a “peak” of the deltoid through the angle at its base. The fantom 
solution z\ = Z3 = 00, z 2 = Z4 is also present in this case. 

3 . 1 . 3 . The isogram. Under assumption a\ = 03 7^ 02 = 04 equation ^ 
becomes 

(ai - a 2 )zfz 2 + 2a 2 z\z 2 - (04 + a 2 ) = 0 
which factorizes as 


(ziz 2 - l)((ai - a 2 )z!Z 2 + (ai + a 2 )) = 0 

The configuration space consists of two lines: z 2 = zf 1 , consisting of paral¬ 
lelograms, and z 2 = zf 1 , consisting of antiparallelograms. 


3 . 2 . Conic quadrilaterals: parametrization by trigonometric func¬ 
tions. If equation ( 12 ) has exactly one solution, then two cases must be 
distinguished: either the two pairs of opposite sides add up to the same 
total length, or two pairs of adjacent sides do so. 

In what follows we will be repeatedly using the fact that a + b = c + d 
implies 

ac — bd = (a + b)(a — d), cd — ab = (a — c)(a — d) 


3 . 2 . 1 . Circumsribable quadrilaterals: 01 + 03 = a 2 + 0,4. 


Proposition 3 . 2 . Assume that 01 + 03 = a 2 + 04 is the unique solution of 
the equation © and that 03 = o m i n . Then the configuration space Qg r (a) 
is isomorphic to the one-point compactification of C/ 27 rZ, and a bijection 
C/ 27 tZ Q° r (a) \ {00} can be established by the parametrization 


zi=pisinf, z 2 — p 2 sin(f -(- i<j ), 


7r 


z 3 = P3 sin [t — — , Z4 = P4 sin [t + — + ia 


n 


where 


Pi = 


(ai + 03) (ai — a 2 ) Faiaf 


02 03 


- 1 


0203 


with p 2 , p^, P4 obtained by cyclically permuting the indices, and 


a = In 


0-2 ®4 


0204 — 0103 


0103 


O2O4 — Ol ®3 


The square roots of negative numbers are assumed to take value in il 


Proof. Due to ai + 03 = <22 + 04 equations Q and ([8]) become 
( 13 a) — 02032^ + aia4z| + (0404 — 0203) = 0 


( 13 b) 


02(01 - 04)2:4 + a 4 (ai - a 2 )zl 


2a 2 a4Z\z 2 + 04(04 + 03) = 0 


Equation ( 13 a) rewrites as 
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with 


Pi = 


a 104 
02^3 


- 1 , 


P 3 = 


Q2«3 

0404 


-1 


so that the coordinates (21,23) can be parametrized as 

21 = p\ sinf, 23 = p 3 cos t 


The substitution of 21 = p\U\ and 22 = p 2 u 2 in ( 13 b) under the assump¬ 
tion ai = a max results in 


u\ + u\ — 2 


0204 


-U1U2 + 


0103 


= 0 


O2O4 — O1O3 0204 — O1O3 

which has the form 

( 14 ) u 2 + v 2 — 2 cos ritu — sin 2 r = 0 

with r = ia and a as stated in the theorem. On the other hand, u and v in 
( 14 ) can be parametrized as u(t) = sinf and v(t) = sin(f + r). This leads to 

21 = pi sin t, 22 = P 2 sin(t + ia) 

The relation between each pair of opposite angles imply that the above 
parametrization is completed by 


Z 3 = P 3 sin (t ± - 


7r 


24 = p 4 sin ( t + ia ± — 


7 T 


for some choices of the signs. To determine the signs, one may analyze 
the equations between 22 and 23 and between 23 and 24 in a similar way. 
Alternatively, one looks at special configurations of quadrilaterals. 

□ 


3 . 2 . 2 . Edge lengths satisfy 01 + 02 = 03 + 04. Equation ([T]) relating 21 and 
23 takes the same form as in the case of curcumscribable quadrilaterals: 

—a 2 o 3 2 2 + 01042! + (oia 4 - a 2 a 3 ) = 0 

so that we have again a parametrization 

21 = pi sint, 23 = p 3 cos t 

But for the pair ip 2 , +4 we have 

(0304 — aia 2 )z 2 z 2 — a\a 2 z 2 + 03042! = 0 

or, equivalently, 

— 03042^ 2 + aia 2 zf 2 + (0102 — 0304) = 0 

which leads to 

zf 1 = p 2 sin t, zf 1 = P4 cos t 

Similarly to the preceding section, we obtain the following (note a sign 
change for 23). 
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Proposition 3 . 3 . Assume that ai + <22 = a 3 + 04 is the unique solution of 
the equation © and that a 3 = a m i n . Then the configuration space Q^ r (a) 
is isomorphic to a one-point, compactification of C/ 2 irZ, and the bijection 
C/ 27 tZ —¥ Qor( a ) \ {°°} can be established by the parametrization 


z 1 = pi sin t, z 2 1 = P2 sin(t + icr ), 


z 3 = p 3 sin (t + - 


7 T 


-1 


7 r 


= p4 sin (t + — + icr 


where 


Pi = 


Oi(l4 

a.2 «3 


- 1 


with p2, p3, pa obtained by cyclically permuting the indices, and 


a = In 


Cl2®4 


<220-4 — <2-1 <3,3 


+ 


0103 


(l2®4 — 


we 


3 . 2 . 3 . The real part of the configuration space. In Theorems 3.2 and 3.3 
have pi,P2 £ M and P3,P4 G iM. It follows that the variables Zi take real 
values if and only if Re t = f (mod7r). The parameter change t = % + iu 
leads to a parametrization that makes the real part of Q° r more tangible. 

Corollary 3 . 4 . Assume that <21 + 03 = <22 + 04 is the unique solution of the 
equation ( 12 ) and that a 3 = a m ; n . Then the configuration space Q° r can be 
parametrized as 


( 15 ) 


where 


and 


z\ = pi cosh u 
Z2 = P2 cosh(it + a) 

Z3 = ip3 sinh u 

Z4 = —ipi sinh(it + a) 


Pi = 


<2404 

a.20.3 


— 1 etc. 


cosh a = 


«2®4 


<220-4 — 0403 


sinh cr = 


0403 


a 2«4 — 0403 


See Figure [ 5 | 

If ai + a 2 
there is a parametrization 


If ai + a2 = 03 + 04 is the unique solution of (12), and a 3 = a m i n , then 


( 16 ) 


z\ = pi cosh u 
zf 1 = P2 cosh (it + a) 
z 3 = —ip 3 sinh u 

zf 1 = —ip4 sinh(u + a) 


with pi and a computed by the same formulas. See Figure [|| 


3 . 3 . Elliptic quadrilaterals: parametrization by elliptic functions. 
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Figure 5. The real part of the configuration space in the 
case cl\ T 113 = a 2 T (14. 



Figure 6 . The real part of the configuration space in the 
case cl\ T C 12 ~~ a 3 T u 4 . 


3.3.1. Special biquadratic equations and elliptic curves. A biquadratic poly¬ 
nomial is a polynomial of degree 2 in each of its two variables. Here we will 
look at a special class of biquadratic polynomials. 

Lemma 3.5. Let m € C \ {0,1}. Then the equation 

(17) u 2 + v 2 = 1 T mu 2 v 2 

defines an elliptic curve. 


Proof. The substitution 


x = u, y = v(l — mu 2 ) u = x, v = 


V 


1 — mx 2 


transforms equation into 

y 2 = (1 — x 2 ){l — mx 2 ) 


If m {0,1}, then the polynomial on the right hand side has four distinct 
roots, thus the equation defines an elliptic curve. □ 
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It is possible to transform equation 0 to the form © by a substitution 
z\ = piui, Z3 = P3U3. In order to make the formulas more compact we need 
some handy notations. 

Let a, b,c,d£ M. Put 


( 18 ) 


s := 


a b c ~\~ d 


a := s — a = 


—a + b + c + d 


b := s — b, c := s — c, d := s — d 

We have a + b = c + d, s — a — b = — (s — c — d) etc. Other, less obvious 
identities, are collected in the following lemma. 


Lemma 3.6. In the above notations we have 

ab — cd = (s — a — c)(s — b — c) 
ab — ab = s(s — c — d) 

abed — abed = s(s — a — b)(s — b — c)(s — a — c) 

Now equation Q can be rewritten in any of the following two ways: 
(19a) (s — a 3 — 04 ) (s — 02 — a^zfz^ — did±z\ + 0203 ,?! + s(s — 02 — 03 ) = 0 

(19b) (0203 — dia,4)zfz^ — didizf + 0203^3 + (0404 — 0104 ) = 0 


Proposition 3 . 7 . The substitution z\ = piu, Z3 = P3V with 

( 20 ) 


Pi = 


's(s - a 2 - a 3 ) 


0404 


0404 

0404 


- 1 , 


P 3 = 


' s(s — a 1 — 04) 

0203 


02^3 

a 2 03 


- 1 


= 1 - 


01020304 


transforms equation 0 to © with 

( 21 ) m 02 03) (s 02 04) (s 03 04) 

01020304 04020304 

In particular, if ( 12 ) has no solutions, then Q^ r (a) is an elliptic curve. 

Proof. An equation of the form 

Az 2 w 2 + Bz 2 + Cw 2 + D = 0 

with D / 0 can be brought into the form © by the substitution 

2 LI 2 2 ^2 

* = --g« , W = -—v 


This yields m = Applying this to the equation ( 19 a) and taking into 
account the third identity of Lemma 3.6 gives us the value of m as stated in 
the proposition. Note that B, C / 0 due to ( 4 b), and D / 0 , since otherwise 

04 + 04 = 02 + 03. 

Further, rn 1 because a* / 0 , and m /0 because this would result in 
s = ai + aj for some i, j. □ 
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3.3.2. Biquadratic equation and Jacobi elliptic functions. Jacobi’s elliptic 
functions can be defined for any complex value k ^ {0,±1} of the Jacobi 

modulus, [ 21 , 33 - 

Lemma 3.8. Let m E C \ {0,1}. Then the curve 
(22) x 2 + y 2 = mx 2 y 2 + 1 

possesses the parametrization 

x(t) = sn(t; k ), y(t ) = sn(f + K]k) with k = \fm 

as well as the parametrization 

x(t) = cn(f; k), y ( t ) = cn(t + K; k) with k = 



Proof. The sn-parametrization follows from the identity 




1 — sn 2 (t;k) 

1 — k 2 sn 2 (f; k) 


The cn-parametrization follows from the identity 


cn(t + K ; k ) 


z sn(t; fc) = , f 1 — cn 2 (f; fc) 

dn(t; &) ' y (A ;') 2 + k 2 cn 2 (t; fc) 


□ 


Remark 3 . 9 . By a parameter change or by symmetry reasons one can 
obtain alternative parametrizations of ( 22 ): 

x = snt, y = sn(t — K) and x = cnf, y = cn(t — K) 


The value of m in our case is given by (21 ), so that m E (—oo, 0 ) U ( 0 , 1 ). 
Let us make a case distinction in order to bring the Jacobi parameter k into 
its usual range ( 0 , 1 ). 


Proposition 3 . 10 . Let a m i n and a max denote the minimum, respectively the 
maximum value taken by 01,02,03,04, and let a m \n + a max s. 

( 1 ) If a m i n + a max < s, then the complexified configuration space Q^ r (a) 
is isomorphic to the quotient of C by a rectangular lattice, which is 
the period lattice of the elliptic sine function with the parameter 

04020304 
04020304 

The tangents of halves of the opposite angles <pi and 933 can be 
parametrized as 

Zi = pi sn(t; k), z 3 = p 3 sn(t + K ; k) 
with pi and p 3 as in (p0|) . 


( 0 , 1 ) 
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(2) If a m m + fl im y > s, then the complexified configuration space Q^ r (a) 
is isomorphic to the quotient of C by a rhombic lattice, which is the 
period lattice of the elliptic cosine function with the parameter 


k = 


t _ 01020304 £ ((J1) 

O1O2O3O4 


The tangents of halves of the opposite angles p i and ip 3 can be 
parametrized as 

zi = pi cn(t; k), z 3 = p 3 cn(t + I \; k) 
with pi and p 3 as in 


Proof. This is a direct consequence of Proposition 3.7 and Lemma 3.8 The 
only things to observe are that a m i n + a max s if and only if ( 12 ) has no 
solutions, and 


Omin + ®max < S s(s - Oi - 0 2 )(s ~ a 2 ~ a 3 )(s - a 3 - Ol) < 0 

44 - 01020304 < 01020304 


□ 


3.3.3. Shift parametrizations of more general biquadratic equations. A gen¬ 
eral symmetric biquadratic equation can be parametrized by a shift of an 
elliptic function: x = f(t), y = /(t + r). This is essentially due to Euler 1 11 
who used biquadratic equations to integrate the relation 


dx dy 

7W) + 7W) 


deg P = 4 


Euler’s argument in a special case is reproduced in [23]. The following lemma 
is a special case that we need in order to parametrize equation ([ 8 ]). 


Lemma 3.11. The curve 

(23) x 2 + y 2 = Ax 
possesses the parametrization 

x(t) = sn(t; k), 

if k £ (0,1) and t£C satisfy 

(24) A = k 2 sn 2 r, B = 
and the parametrization 

x(t) = cn (t), 

if k E (0,1) and r E C satisfy 

(25) A = -k B = 

dn r 


2„,2 


y + 2 Bxy + C 


y(t) = sn (t + t; k) 


cn r dn r, C = sn 2 t 


y{t) = cn (t + t) 


dn 2 t 


, C = {k'f 


dn 2 r 


Here k! = y/l — k 2 is the conjugate modulus. 
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sntcnr dn r + cn t dn t sn r 
1 — k 2 sn 2 t sn 2 r 


Proof. The addition law 

sn (t + t) = 
can be rewritten as 

(1 — k 2 sn 2 t sn 2 r) sn(t + r) — cnr dn r sn t = \/(l — sn 2 t)( 1 — k 2 sn 2 f) 

By taking the squares of both parts and substituting x = snt, y = sn(f + r) 
we obtain 

k A sn 4 r x 4 y 2 — k 2 sn 2 r x 4 — 2A ; 2 sn 2 r x 2 y 2 + 2 k 2 sn 2 r cn r dn r x 3 y 

+ (1 + k 2 sn 4 r)x 2 — 2 cn t dn t xy + y 2 — sn 2 t = 0 

The same equation holds with x and y exchanged. After antisymmetrizing 
and dividing by A; 2 sn 2 r(x 2 — y 2 ) we obtain the equation (23). 

In the cn case the proof is similar. □ 


A biquadratic equation (23) can be viewed as an implicit form of addition 
formulas for sn or cn. 

3.3.4. Parametrization theorems. We are ready to derive a simultaneous 
parametrization of all four angles of a quadrilateral with fixed side lengths. 
With the help of notation (18) equation Q can be rewritten as 

(26) ai(s — a\ — a 3 )zjz 2 + 04(5 - a 3 - a 4 )zf + a 2 (s - a 2 — a 3 )z 2 

— 2a 2 a 4 ziz 2 + 503 = 0 

Cyclic shifts of indices produce equations for other pairs of adjacent angles. 

Proposition 3.12. Let a = (ai, a 2 , a 3 , a 4 ) be such that ai±a 2 ±a 3 ±a 4 ^ 0 
for all choices of the signs. Define 


Pi = 


's{s - a 2 - a 3 ) 
0104 


Gq <24 

aia 4 


- 1 £ 


U*I 


and similarly p 2 , p 3 , and p 4 through a cyclic shift of indices. The configura¬ 
tion space Qg r (a) has one of the following parametrizations in terms of the 
tangents of the half-angles of the quadrilateral. 

(1) If cimin T Umax ^ s and U 3 = u m ; n , then 
z 1 = Pi sn t, z 2 = p 2 sn (t + r), z 3 =p 3 sn (t - K ), z 4 = p 4 sn (t + K + r) 

where the modulus of sn is 


k= /1 _ ^304 6{01)i 

V ^0^304 

and the phase shift r is determined by 

dn t = \ I a2a _ 4 , r £ (0, iK') 
a 2 a 4 
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(2) If ®min “I - ®max ^ S and d\ — Umax; then 
Zi=picnt, z 2 = p-2 cn(t + t), z 3 = p 3 cn(t - K), Z 4 = p 4 cn(t + K + r) 
where the modulus of cn is 


k = 


!- ° 1<Wl4 6 ( 0 ,!), 


(21020304 

and the phase shift r is determined by 


dnT = 


O 20 4 

02 o 4 


r G (0, iK') 


Proof. Substituting zi = and Z 2 = P 2222 and using the identities from 


Lemma 3.6[ we transform (|26|) to 
(27) u\+U 2 = 


— 2 


020 4 


0102 -- 


0103 


020 4 * \/a2a 4 \/a20 4 — 0103 020 4 — 0303 

The range of dn on (0 ,iK') is (l,+oo). Due to Lemma 


3.6 


we have 


020 4 > 020 4 , thus there exists r G (0 ,iK') with dnr = ■ From this 

we compute 


sn t = — 


a\a 3 


a 2a 4 — 0103 ’ 


cn t = 


0 , 20,4 


02O4 — 0403 


Besides, cnr > 0 for r G (0, iK'). It follows that the coefficients of (27) 


coincide with the coefficients of (23) for given values of k and r. A similar 


check can be made for other pairs of adjacent angles (where one should be 
careful with extraction of square roots of negative numbers), which finishes 
the proof of the sn-parametrization. 

□ 


In the case a m i n + a max > s use the second part of Lemma 3.11 


3.3.5. The real part of the configuration space. Under the assumptions of 


Theorem 3.12 we have p \, p 2 G R-|_, p 3 ,P 4 G It follows that all of the 

variables Zi take real values if and only if 

Ret = 7i(mod2/\), for the sn-parametrization 
Ret = 0(mod2/\), for the cn-parametrization 

As in the trigonometric case, we make a variable change so that real values 
of parameter yield real angles pi. 

In the sn-case put t = iu+K , and in the cn-case put t = iu. Jacobi’s imag¬ 
inary transformations [31 Chapter 2] allow to rewrite the parametrizations 


of Theorem 3.12 in terms of elliptic functions with the conjugate modulus. 
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Proposition 3.13. (1) Let a m i n +a max < s, and assume that a 3 = a m j n . 

Then the configuration space Q° r can be parametrized as 


z 1 = 


Pi 


Z3 = 


dn(it; k') ’ 
P3 


Z2 = 


P2 


. sn(it; k') 
dn(tt + iK ; k') P3 cn(u; k' 


z 4 = 


Pa 


dn(rt + cr; k') 

■n(w; k') 

:n(rt; k' 
sn (u + a; k') 


dn(u — iK + cr; k') 


= ~IPA 


cn {u + a; k') 


where 


k , = (01020304 (0 ,l )j 

Q.lfl2®3®4 


the amplitudes pi are given by 


Pi = 


! s(s - a 2 - a 3 ) 


a.1 <2.4 
aia4 


— 1 , etc. 


and the phase shift a is determined by 
sn(cr; k') = 


0403 


, (T€{0 ,K') 

a 2 a4 

(2) Let a m i n + a max > s, and assume that a\ = a mE 
ration space Q^ r can be parametrized as 

Pi P2 


Then the configu- 


zi = 


z 3 = 


z 4 = 


cn (it; k') ’ 

P3 

cn(u + iK: k') 
Pa 


Z2 = 


cn(it + a; k') 
_ , sn 0;fc') 

' P3 dn {u-k') 


cn(u — iK + a) 


sn(u + cr; k') 

= -ikp 4 —— - 777 

dn(tt + cr; k') 


where 


k' = 


04020304 


04020304 
the amplitudes pi are given by 


e (0,1), 


pi = 


1 s(s - 02 - 03 ) 
0404 


0404 

0404 


— 1, etc. 


and the phase shift 0 is determined by 


sn(ir; k') = / QlQ3 ? u e (0,K') 
V 0204 


The shapes of quadrilaterals corresponding to special values of parameter 
t are shown on Figures [7] and [8} Note that the real part consists of one 
component (a quadrilateral can be transformed to its mirror image by a 
continuous deformation) if a m i n + a max > s, and of two components (no 
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Figure 7. The real part of the configuration space in the 
Case Omin T ®max *7 S. 



Figure 8 . The real part of the configuration space in the 
Case Umin T a max > S. 


continuous deformation between a quadrilateral and its mirror image) if 

®min T flmax ^ S. 

Since = cot < 7‘. Proposition 


3.13 


implies Theorem pi 


4. The configuration spaces of non-oriented quadrilaterals 


4.1. Double covers. For a = (01,02,03,04) satisfying the conditions (4a) 
and (4b), let Q(a) denote the space of congruence classes of quadrilaterals in 
M 2 with side lengths a in this cyclic order. This time, orientation-reversing 
congruences are allowed, so that to almost every element of Q(a) there cor¬ 
respond two elements of Q or (oi). Exceptions are quadrilaterals with collinear 
vertices. 

If a is elliptic, that is ai ± 02 ± 03 ± 04 / 0 for every choice of signs, then 
there are no exceptional quadrilaterals and the natural map Q or (a) —> Q(a) 
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is a double cover (disconnencted for non-Grashof quadrilaterals). The equiv¬ 
alence relation on Qor(«) looks very simple in terms of the angles of a quadri¬ 
lateral: 

(<Pi, <P2, <P3, Pa) ~ (~Pi, ~P2, ~P3, -‘Pa) 
and has the same form in the variables z% = tan j . By extending this 
to Q& r (a) we dehne the complexified configuration space of non-oriented 
quadrilaterals 

Q C (a) := Q^.(a)/(zi,z 2 ,z 3 ,z 4 ) ~ (~z 1: -z 2 ,-z 3 ,-z 4 ) 


The quotient map 

Qor( a ) Q C ( a ) 

is a double cover between elliptic curves. In terms of the holomorphic pa¬ 
rameter t of Proposition 3.12 we have t ~ t + 2 K both in the sn and in the 
cn case. This identifies Q{a) with C/2/iZ + 2K'iZ and proves the first part 
of Theorem [3l 


4.2. Diagonal lengths as coordinates on the configuration space. 

The cosines cos cpi are meromorphic functions on the elliptic curve Q c (a). 
On the other hand, by the cosine law they are linearly related to the squares 
of diagonal lengths. Since a quadrilateral with fixed edge lengths is uniquely 
determined up to congruence by the lengths of its diagonals, this provides a 
natural embedding of Q c (a) in (CP 1 ) 2 . 

Lemma 4 . 1 . Let x and y be the lengths of diagonals (x is separating a\ and 
a 2 from a 3 and a 4 ), and 

2 2 
u = x , v = y 

Then the configuration space Q(a ) is an algebraic curve given by the equation 
(28) u 2 v + uv 2 + 2d\\uv + d\ou + d^v + doo = 0, where 

d\\ = —-(«i T T a 2 + a 2 ) 

dio = (a? - a 2 )(a| - a\) 

doi = (a? - a|)( a 4 - °i) 

doo = (af - a\ + a§ - a\)(a\al - a 2 2 a\) 


Proof. Viewing a planar quadrilateral as a tetrahedron of zero volume, we 
equate to zero the Gram determinant of the vectors (two sides and a diago¬ 
nal) issued from one of its vertices. 


det 


^ a i 

a^-\-x 2 — a| 

i af+ai-y 2 
\ 2 


a 2 -\-x 2 — a 2 
2 

x 2 


x 2 +a 2 — a 2 
2 


a i+ a I -y 2 ' 

x 2 +af— a| 


a 4 


= 0 


A simple computation leads to the equation in the lemma. 


□ 
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4.3. Conjugate quadrilaterals. Recall the following notations. 
Definition 4 . 2 . For a € M 4 the conjugate quadruple a 6 M 4 is defined as 


—a i + a2 + 0.3 + 04 
ai := s - ai =--- 

ai + 02 — 03 + U4 
«3 : = s — 03 =---, 


fli — fl2 4 (13 4 ^4 
02 := s - 02 = --- 

_ a\ + 02 4 “ 0-3 — a 4 

04 := s — a 4 = -—- 


Lemma 4 . 3 . The map a ^ a is an involution on the space of quadruples 
a G M 4 satisfying inequalities pal) and (4b). 


Proof. We have a = a since ^2 at = 22 , a i- The quadrangle inequality s > a* 
is equivalent to the positivity of the length a* > 0 . □ 


Thus a quadruple conjugate to the edge lengths of a quadrilateral can 
itself be used as edge lengths of another quadrilateral. 

Note that a = a if and only if all ai are equal, and a coincides with a up 
to the dihedral group action on { 1 , 2 ,3,4} if and only if a m ; n + a max = s. 
Thus for ai ± 02 ± 03 ± 04 / 0 quadrilaterals in Q{a) are not congruent to 
quadrilaterals in Q(a), even if we ignore the marking of sides. 

The more surprising is the second part of Theorem [3j for every quadrilat¬ 
eral with side lengths a and diagonal lengths x and y there is a quadrilateral 
with side lengths a and diagonal lengths x and y. For this, it suffices to 
show that the coefficients dy in (28) don’t change if we replace a by o. This 
follows from a + b = c + d, a — b = b — a, which implies that the coefficients 
doi and dio don’t change and from Lemma |4.4| below, which takes care of 
dn and doo- 


Lemma 4.4. Let 

a + b + c + d _ 

s = -, a = s — a, 

2 

Then we have 

a 2 + & 2 + c 2 + d 2 = -2 + 5 2 + -2 + (? 
ab + cd = ab + cd 
ab — cd = ^(c 2 + d 2 — a 2 — b 2 ) 
ft is easy to see that 


&min + ^max ^ S &min ~\~ fl ma x S — S 

Thus Qo r (a) is a quotient of C by a rectangular lattice if and only if Q° r (a) is 
a quotient of C by a rhombic lattice. Thus the orientation forgetting covers 
over Q c (a) = Q c (a) look as shown on Figure [ 3 J and Theorem [ 3 ] is proved. 









DEFORMATION OF QUADRILATERALS AND ADDITION ON ELLIPTIC CURVES 25 


5. Periodic quadrilaterals 


5.1. Folding of conic quadrilaterals. 


Proposition 5 . 1 . For a quadrilateral of conic type, the composition of fold¬ 
ings 

F 3 oF 4 : QS-( o) _> QC( o) 

has a unique fixed point, namely the “flattened” quadrilateral with collinear 
vertices. 

For any initial position q E Q f f r (a), the sequence q n = (F 3 o Ffl) n (q) 
converges to the flattened quadrilateral, with the angles of q n tending to 0 or 
it exponentially. 


Proof. Proposition 3.2 implies that in terms of the parameter t the folding 
maps take the form 


F 3 (t) = 7T — t, Fflt) = 7r — 2 ia — t 
so that (F 3 o Ffljft) = t + 2ia. Thus 

(F 3 o F 4 ) n (t) -A oo 


and oo is a unique fixed point of F 40 F 3 . The exponential decay of the angles 
follows from 


zi((F 3 o Ffl) n (t)) = pi sin(i + 2 nia) ~ e 2ncr 
and z\ = tan ff ~ —-— as czn —> it. □ 

5.2. Folding of elliptic quadrilaterals. 


Proof of the Darboux porism. Proposition |3.12| implies that the folding maps 
Fi: Q^ r (a) —> Q^ r (a) are involutions of the form t i-A- U — t. Thus the 
composition F 3 o F 4 is a translation t >-A t + (t 3 — t^). The n-th iteration 
(F 3 oF 4 ) n is a translation by n{t 3 — tfl). If it has a fixed point, then n(t 3 — ti) 
belongs to the lattice of periods, and consequently (F 3 o Ffl) n = id. □ 


The definition of n-periodicity given in the Introduction depends on the 
choice of an ordered pair of adjacent vertices. Let us clarify this dependence. 

Lemma 5 . 2 . If the map F 3 o F 4 : Q^ r (a) —> Q^ r (a) has order n, then the 
maps Fi o F 3 , F\ o F-2, and F2 o F\ also have order n. The map F2 o F 3 also 
has a finite order, but possibly different from n. 


Proof. The maps F4 o F 3 and F 3 o F4 have the same order because they are 
conjugate to each other. Proposition 3. 12| implies that the maps F 3 o F4 and 
F\ o F 2 are translations by the same vector (modulo the lattice of periods), 
therefore have the same order. Finally, the same theorem shows that 

F 3 o F 4 (t) = t + 2r, F2 o F 3 (t) = t — 2K — 2 t 


Thus one translation is commensurable with the lattice of periods if and 
only if the other is. □ 
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Definition 5 . 3 . A quadrilateral is called n-periodic, if both F;> o F± and 
+2 o F3 have order n. It is called (n, m)-periodic, if one of the maps F3 o F4 
and F2 o F3 has order n, and the other has order m. 


On the other hand, the folding period modulo orientation , that is the 
order of the natural map F30 F4: Q c (a ) —>• Q c (a ) does not depend on the 
choice of a pair of adjacent vertices because F\ = F 3 and F2 = Fa- 


Proposition 5 . 4 . Let a E ( 0 ,/i 7 ) be the basic shift as defined in Theorem 


3 . 13 . Assume that a = —K' with m and n coprime. Then all quadrilaterals 
with side lengths (01,02,03,04) are n-periodic up to symmetry. Besides, 

(1) in the sn -case 


n is even =>■ the period length is n 
n is odd =>■ the period lengths are (n, 2 n) 


(2) in the cn -case 

n is even =>■ the period length is 2n 
n is odd =>■ the period lengths are (n, 2ra) 


Proof. Let us use the parameter change made in Section [3.3.5 see Figures 
[7] and [8] for illustration. In terms of the real parameter u the two shifts 
performed by different pairs of adjacent foldings are 


2(j and 2 iK — 2 a 


The n-fold iteration results in the shifts by 

n ■ 2 a = 2 mK', n( 2 iK — 2 a) = 2 niK — 2 mK' 

These are the smallest multiples of 2<r that belong to the lattice 2/i , Z+2iA"Z, 
which is the period lattice for Q c (a), see Section 

More details are needed to determine the period on Q^ r (a). If n is even, 
then m is odd. Thus both of the above shifts are periods of dn and half¬ 
periods of cn. If n is odd, then the first shift is a period of dn, but the second 
only a half-period. If m is also odd, then the first one is a half-period for 
cn, and the second a period; vice versa if m is even. □ 


4.1 


5.3. Quadrilaterals with small periods. 

5.3.1. Period 2 : orthodiagonal quadrilaterals. The composition of foldings 
F 3 o F 4 has order 2 if and only if F 3 and F 4 commute. This happens if and 
only if the third and the fourth vertices stay on their respective diagonals 
after the folding. This, in turn, is equivalent to diagonals being orthogonal. 
Pythagorean theorem implies that diagonals of a quadrilateral are orthog¬ 
onal if and only if 0^ + 03 = 02 + 04. This proves the period 2 case of the 
Darboux porism without using elliptic functions. 
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B 



Figure 9. Folding of an orthodiagonal quadrilateral. 

Proposition 5.5. A quadrilateral is 2-periodic if and only if its diagonal 
are orthogonal to each other. This property depends only on the side lengths 
of the quadrilateral, and is equivalent to 

(29) 


a ^ “l - U 3 — a<2 d - dq 


Note that in the orthodiagonal case equation ([ 8 ]) can be rewritten as 
{AlZl + B\Z 1 !)(A2Z2 + B 2 Z 2 1 ) = 1 
which also makes the commutation of foldings apparent. 

5.3.2. Period 4. Introduce the notation 

A = Q. 1 U 3 , B = A = aid 3 , B = a 2 O 4 


Then by Theorem 3.13 we have 


AB 2 A 


(30a) (k') 2 = -=~^, sn 2 7 = ^ 

(30b) (k') 2 = 


AB' 

AB 

AB' 


B ’ 

2 A 
sn a = - , 
B 


9 , 9 

cn (7=1 — —, dn (7 = 1 

JD 


2 1 A 

cn (7=1-- , 

B 


dn 2 7=1- 


A 

B 

A 

B 


in the sn- and in the cn-case, respectively. (All elliptic functions have elliptic 
parameter k '.) Note that the numbers A,B,A,B satisfy the identities 

A + B = A + B, AA- BB = (A + B)(A- B) 

Also we have a m i n + a max < s AB < AB. 


From (30a) we compute 


sn 2 u = 


2sn7cn7dn7 2 V. AB 


1 — k 2 sn 4 7 


A + B 


and only if sn(27; k') = 1, which is equivalent to 

2 \[AB = A + B 


This gives another proof of Proposition 5.5 since a £ (0, K') equals if 


that is, to A = B, which by Lemma 4.4 is equivalent to a 2 + + a z A . 
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Proposition 5 . 6 . A quadrilateral is 4 -periodic if and only if one of the 
following takes place: 

( 1 ) flmin ~t~ ®max ^ S and CliO3 = O2O4/ 

( 2 ) fl m i n T Umax ^ S and 

(0103 — 0204) (0103 + 0204) = ±2^/04020304(04 + a 2 — a 2 — a 2 ) 


Proof. According to Proposition | 5 . 4 t a quadrilateral is 4 -periodic if either 
Omin + Omax < s and a E {%-, or o min + a n 


> s and a = r f 2 . 


In the latter case by symmetry with the sn-case we have 


a = 


K' 2 y/AB 




2 "A + B 
In the former case we compute 


= 1 ^ A = B 


sn2cr = 


2VAB 

A + B 


cn 2 a = + 


\A-B\ 
A + B 


dn 2a = 


\A~B\ 
A + B 


where cn 2 o > 0 if and only if a E (0, Using once again the formula for 
sn2x, we obtain 


sn 4 <r = ±4 


y/AB\A-B\ 
(A + B)\A - B\ 


Since a E {^p, ^-} sn 4 a = +1 , we obtain the second condition. 


□ 


5 . 3 . 3 . Period 3 . 


Proposition 5.7. A quadrilateral is ( 3 ,6)-periodic if and only if 

(0103 + 0204) 2 E {404030204,402040103,401030403,402040204} 

Proof. Due to Proposition | 5 . 4 | we have to find the necessary and sufficient 
conditions for a E -}. For a E ( 0 ,K') we have 

K' . . . cn a 

a = — sn2cr = sn(A — a) = - - 

3 dno 

2 K' 

a = —— 4 => sn 2 ct = sn(2 K — a) = sn a 

This leads to the equations (A + B) 2 = 4 AB and (A + B) 2 = 4 BB in the sn- 
case, and to (A + B) 2 = 4 AB and (A + B) 2 = 4 AA in the cn-case. Since we 
were making assumptions 03 = a m i n and ai = a max in the sn- and cn-cases 
respectively, omitting them leads to symmetrizing with respect to A and B. 
Thus all four conditions are valid both in the sn- and in the cn-case. □ 

Example 5.8. For 01 = 02 = 03 = a and 04 = b we have 04 = 02 = 03 = 
and 04 = so that 

(0103 + 0204) 2 = a 2 (a + b) 2 = 401030403 

and hence every quadrilateral with side lengths (a, a, a, b) is ( 3 , 6)-periodic. 
For a geometric proof of this fact, see m Figure 7 ]. Note that the quadri¬ 
lateral is Grashof for a > b and non-Grashof for a < b. 
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5 . 4 . The general periodicity condition. We have a = ffK' (and thus 
the quadrilateral is periodic) if and only if 

cn(ncr; k ') E {0, — 1,1} 

It is possible to express cn(ncr) as a rational function of cncr, so that 
<7 = ^ K' becomes an algebraic conditions on the side lengths. These ra¬ 
tional functions can be computed recursively, with the help of the following 
formulas that generalize the Chebyshev polynomials. 

/?,,(/) = 1 , Rl(t)=t, Rn+l(t) = 1 _ - Rn-l(t) 

Benoist and Hulin [ 5 ] relate the dynamics of quadrilateral folding to the 
Poncelet dynamics on the circles of radii R = 0204 and r = ^(a|— a^+a^— of) 
whose centers are at the distance a = 0403 from each other. From the known 
conditions on short trajectories they obtain characterization of quadrilater¬ 
als with small folding period. One may then use Cayley’s explicit conditions 
on the periodicity of Poncelet trajectories to obtain the periodicity condi¬ 
tion for quadrilateral folding. We choose a more direct approach, following 
the general principle derived by Griffiths and Harris m that also leads to a 
modern proof of Cayley’s conditions. This results in our Theorem [ 4 ] which 
we prove below. 


Theorem 5 (Griffiths-Harris [ 16 ]). Let E be an elliptic curve andp,q £ E. 
For an integer n > 3 we have n(p — q) = 0 with respect to the additive 
structure on E if and only if the following condition is satisfied. 

Choose a meromorphic function x of order 2 on E such that x has a 
double pole at p and a zero at q. Let a, b, c £ C be the values of x at its 
branch points other than p, so that E is isomorphic to the Riemann surface 
of the function 


y 

'e' 

1 

> 

II 

on p 

- q being of 

A2 

A 3 

^3 

A4 

Ak+i 

Ak+2 ■ ■ • 

A3 

a 4 

A\ 

a 5 

Ak+ 1 

Ak+2 ■ ■ • 


^ 4 fc+i 

Ak+2 

A-2 k 

Ak+i 

Ak+2 

A2k-1 


i =0 


= 0, if n = 2 k + 1 


= 0, if n = 2 k 


Our first goal is to construct a function x as described in the Griffiths- 
Harris theorem. Consider the function u equal to the squared length of the 
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Figure 10. The quadrilaterals p and q = F% o F4(p). 


diagonal separating the sides 04,02 from 03,04. At the branch points, u 
takes the values 

(ai + a2) 2 , (oi 02) 2 , (o3 + a4) 2 , (03 — 04) 2 

Let p € Q{a) be such that u(j>) = (03 + 04) 2 . If 01 + 02 > 03 + 04, then p 
lies in the real part of Q(a) and corresponds to the quadrilateral in which 
the sides 03 and 04 are aligned to form a segment of length 03 + 04. The 
quadrilateral p is a fixed point for the folding F4. and the folding F3 maps 
it to the quadrilateral on Figure [To| right. 

Denote q = F3oF 4 (p). We have (F$ o +4)” = id if and only if n(p — q) = 0 . 
The function 

( 31 ) x= u-ui^ 

u — u{p) 

has a double pole at p and a zero at q, and the same branch points as u. 

Now it remains to compute the values of x at the three remaining branch 
points: 

(oi + a 2 ) 2 - u{q) (01 - a 2 ) 2 - u(q) (a 3 - a 4 ) 2 - u(q) 


(01 + a 2 ) 2 - (a 3 + a 4 ) 2 ’ (04 - a 2 ) 2 - (a 3 + a 4 ) 2 ’ (a 3 - a 4 ) 2 - (a 3 + a 4 ) 2 

Lemma 5 . 9 . We have 

a 3 (a 2 - a 4 ) + a 4 (a 2 - a\) 


v{p) = v{q ) = 

= (o 4 -Q3)(aj ~ a!) + 


03 + 04 

(af — o| + a| — a 2 )(afa| — a 2 a 2 ) 


03 + 04 (03 + 04)(03(a 2 - af) + 04(02 - a|)) 

Proof. By the cosine law we have 

v(p) = a 2 + a 2 + 20404 cos +i(p) 

a\ = a 2 + (a 3 + o 4 ) 2 + 2oi (a 3 + a 4 ) cos (^1 (p) 

Expressing cos <p\ (p) from the second equation and substituting into the first 
gives us the formula for v(p). 


To compute u(q), rewrite equation ( 28 ) as a quadratic equation in u: 


vu 2 + (v 2 + 2dnv + d\o)u + (doiv + d 0 o) = 0 
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Figure 11 . A pair of conjugate quadrilaterals. 

At v = v(p) this equation has two solutions u = u(p ) = (03 + a 4 ) 2 and 
u = u(q). Therefore 

, . \ 2 1 \ d m v(p) + d 00 

(a 3 + a 4 ) u(q) = --- 

v(p) 

and we arrive at the formula for u(q) stated in the lemma. □ 

Finally, a lengthy calculation yields the following formulas for the values 
of x at the branch points: 

(aia 3 + oqgq ) 2 (0103 — 0204)^ (<210,3 — o^cq ) 2 

u(p)v(p) ’ u(p)v(p) ’ u(p)v(p) 

Scaling x by an appropriate functions we obtain a function with the same 
poles and zeros and with the values 

(oiu 3 + gqg^) 2 , (0103 — 020,4)^, (0,40,3 — d 2 U 4 ) 2 

at its branch points. This proves Theorem [ 4 j 

6 . Applications 

6 . 1 . Ivory’s theorem. The second part of Theorem [ 3 ] says that for every 
quadrilateral with the side lengths a and diagonal lengths x and y there is a 
quadrilateral with the side lengths a and the same diagonal lengths x and y. 
We call these quadrilaterals conjugate. See Figure [TT] for an example, where 
we have chosen 

a = ( 10 , 5 , 6 , 3 ), a = ( 2 , 7 , 6, 9 ) 

Arseniy Akopyan has noted [lj that the existence of a conjugate to each 
quadrilateral is equivalent to the Ivory theorem. 

Theorem 6 (Ivory’s Theorem). Take two ellipses and two hyperbolas from a 
confocal family of conics. Then the diagonals of any curvilinear quadrilateral 
cut out by them have equal lengths. 

Proof. Consider a quadrilateral with vertices at the foci and at the end¬ 
points of one of the diagonals of the given curvilinear quadrilateral. The 
order of vertices is chosen so that the foci form a pair of opposite vertices. 
The quadrilateral formed in the same way but using the second diagonal of 
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the curvilinear quadrilateral has side lengths conjugate to those of the first 
quadrilateral. Since this pair of conjugate quadrilaterals has one diagonal 
in common (the segment between the foci), the other diagonals have equal 
lengths. See Figure [l2| □ 




Figure 12. Ivory’s theorem and its proof using conjugate quadrilaterals. 


With the same argument one can prove the Ivory theorem on the sphere 
and in the hyperbolic plane, see Section [7] 

6 . 2 . Dixon’s angle condition in the Burmester mechanism. The link¬ 


age on Figure 13 consists of two small quadrilaterals sharing an edge and 
of one big quadrilateral having an angle in common with each of the small 
ones. It is easy to see that a linkage of this form is rigid in general. The 
following is a necessary condition of its flexibility. 




that is generically rigid, linkage. 


Proposition 6 . 1 . If the linkage on Figure [73] is flexible and both small 
quadrilaterals are of elliptic type, then the angles marked at the left remain 
either equal modulo it or complement each other modulo n during the flex: 

6 \ ± 02 = nir 


Proof. Forget the big quadrilateral and require only that the horizontal sides 
of the small quadrilaterals remain aligned. The configuration space Q 12 of 
this coupling is a fiber product of the configuration spaces Q\ and Q2 of the 
small quadrilaterals, see the diagram (32). 
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(32) Q 12 

Qi Q 2 

IV 1 eCP 1 z e CP 1 CP 1 9 W 2 

Here z = z\ = tan ^ = cot ^ and Wi = tan 

If the branch sets of the maps f\ and $2 on the diagram don’t coincide, 
then Q 12 is a hyperelliptic curve. On the other hand, if the coupling remains 
flexible when restrained by the big quadrilateral, then Q 12 is the configura¬ 
tion space of the big quadrilateral, thus cannot be hyperelliptic. 

Branch points of correspond to 9i E {0,7r}. Thus for flexibility it is 
necessary that 

0] G {0, 7T} <=> 0 2 G {0, 7T} 

Applying Lemma |2.3| to both small quadrilaterals we obtain that cos 9\ and 
cos 62 must be linearly dependent. It follows that either cos0i = cos 62 , that 
is 9\ — 62 = mr, or cos 9\ = — cos 62 , that is 9\ + 02 = mr. □ 

It turns out that if the linkage on Figure [l3| is flexible, then another bar 
can be added without restraining its motion, see [29] . This results in a 
linkage on Figure [14] left, called the Burmester focal mechanism. Of course, 
points on the sides of the big quadrilateral may also lie on the extensions of 
these sides. In particular, if one of them lies at infinity, then the trajectory 
of the central joint is a line orthogonal to the corresponding side. The result 
is a straight-line mechanism of Hart mmm , see Figure [Tfj right. 







Figure 14. Burmester’s linkage and its special case, Hart’s 
straight line mechanism. 

The problem of characterizing all flexible configurations of the form shown 
on Figure [14], left, was studied by Krause [21] using Darboux’s parametriza- 
tion of the configuration space. He found three families: in the first case 
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two opposite small quadrilaterals are parallelograms while the other two are 
similar; in the second all small quadrilaterals are deltoids and the big one is 
orthodiagonal; in the third family all quadrilaterals are conic or elliptic and 
the lengths of 12 segments are subject to 7 polynomial relations. 

More examples of overconstrained flexible linkages can be found in m- 

6.3. Flexible polyhedra. Consider a polyhedral cone made of four faces 
meeting along four edges. If the edges are viewed as ideal hinges, and the 
faces are allowed to cross each other, then the configuration space of the 
cone is that of a spherical quadrilateral obtained by intersecting the cone 
with a sphere centered at the apex. This allows to study the flexibility of 
polyhedral surfaces with vertices of degree 4 by means of elliptic functions. 

The only closed polyhedron with all vertices of degree 4 is the octahe¬ 
dron. Bricard [6] studied the system of biquadratic equations between the 
tangents of dihedral half-angles and described three classes of flexible (self- 
intersecting) octahedra. In fact, it suffices to look at the three dihedral 
angles adjacent to the same face, so Bricard’s method was to find when the 
resultant of two biquadratic polynomials P\(z 2 , 23 ) and P‘ 2 (z \, Z 3 ) is divisible 
by a biquadratic polynomial Ps(zi, 22 ). 

Recently, Gaifullin m classified flexible cross-polytopes (analogs of the 
octahedron in higher dimensions) in the spaces of constant curvature. Unex¬ 
pectedly enough (polyhedra of higher dimensions tend to be “more rigid”), 
flexible cross-polytopes exist in all dimensions, besides, in § 3 there is a 
type that is not present in Bricard’s classification for M 3 . Gaifullin used 
parametrization of dihedral angles by Jacobi elliptic functions, similar to 
the one in Section 13.3.51 




Figure 15. A Kokotsakis polyhedron and the associated 
spherical linkage. 


There are plenty of non-closed 3-dimensional polyhedra with all vertices 
of degree 4. The neighborhood of a face of such a polyhedron consists of an 
n-gon surrounded by a “belt” of triangles and quadrilaterals, and is called 
a Kokotsakis polyhedron , [20]. See Figure 15, left, for a polyhedron with 
a quadrangular base. All polygons are required to remain flat during a 
deformation. Thus the polyhedron is flexible if and only if the spherical 
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right, deforms in such a 


framework with scissors-like joints on Figure 15 
way that the two marked angles remain equal. 

The configuration space of a Kokotsakis polyhedron is described by a 
system of n biquadratic equations in n variables. Kokotsakis EDI, Graf and 
Sauer [ 25 ], and Stachel [ 25 ] described several types of flexible polyhedra. 
In [ 18 | a complete classification in the case of a quadrangular base was 
obtained. Parametrization by elliptic functions is used in an essential way; 
among the new flexible types there is an irreducible one, where the resultant 
of P(z\, Z2) and Q(z2, 23) is an irreducible polynomial and thus proportional 
to the resultant of R(z3, 24) and $(24, 21). 


7 . Configuration spaces of spherical and hyperbolic 

QUADRILATERALS 

Qualitatively, all results of this article apply for quadrilaterals on the 
sphere or on the hyperbolic plane. Configuration spaces of spherical quadri¬ 
laterals play a role in the study of flexible polyhedra, see Section | 6 . 3 | 

For spherical quadrilaterals, the admissible side lengths are subject to 
inequalities 

( 33 a) 0 < ai < n Mi 


( 33 b) 


CLi CLj T CLfe T U/ CLi “t“ 27 T Vi, 


while in the hyperbolic case the conditions are the same ( 4 a), ( 4 b) as in the 
euclidean. Similarly, in place of equation ( 12 ) in the spherical case we have 

ai ± 02 ± 03 ± CI4 = 0 (mod 27 r) 


This adds “antideltoids” and “antiisograms” to the respective cases in Def¬ 
inition 13.11 

In terms of the tangents of half-angles the configuration space of oriented 
quadrilaterals is described by equations 0 and ([8]), where in the spherical 
case we have 

. ai + a,2 — 03 — <24 . «i — CL2 + 03 — «4 

b-22 = sin-sm- 

2 2 

. a\ + 02 + 03 — <24 . 01—02 — 03 — 114 
620 = sm---sm--- 

. fll — 02 + O3 + 04 . Ol + 02 — O3 + O4 

602 = sm---sm--- 

. Ol — 02 — O3 + 04 . Ol + 02 + O3 + O4 

Ooo = sm---sm--- 

and similarly for Cj-j with cn = — sin 02 sin 04. In the hyperbolic case sinh 
takes the place of sin. 

Further, orthodiagonal spherical quadrilaterals are characterized by 


cos 01 cos 03 = cos 02 cos 04 
and these are the only 2-periodic ones. 
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There are analogs of identities from Lemmas 3.6 and 4.4 
are collected in the next lemma. 


Some of them 


Lemma 7 . 1 . For any a, b, c, d the following identities hold. 

cos a cos b + cos c cos d = cos a cos b + cos c cos d 

sin a sin b + sin c sin d = sin a sin b + sin c sin d 

cos a cos b — cos c cos d = sin a sin b — sin c sin d 

sin a sin b sin c sin d + cos a cos b cos c cos d = sin a sin b sin c sin d + cos a cos b cos c cos 

Here a = ^(— a + b + c + d) etc. The same identities hold when sin and cos 
are replaced by sinh and cosh. 

Parametrization of the configuration space in terms of Jacobi elliptic func¬ 
tions (or trigonometric functions, for conic quadrilaterals) are derived in the 
spherical case in J[T8] • 

Equation in terms of diagonal lengths, describing the space of quadri¬ 
laterals up to possibly orientation changing congruences undergoes a more 
significant change. 

Lemma 7 . 2 . Let x and y be the lengths of diagonals in a spherical quadri¬ 
lateral with side lengths 01,02,03,04 (in this cyclic order). Denote 

u = cosx, v = cos y 

Then the configuration space of quadrilaterals up to congruence is an alge¬ 
braic curve given by the equation 

u 2 v 2 — u 2 — v 2 + 2 dnuv + diox + doiy + do 0 = 0, where 
dn = —(cos Ol cos 03 + cos 02 cos 04) 
dio = 2(cos Ol cos 02 + cos 03 cos 04) 
doi = 2(cos oi cos 04 + cos 02 cos 03) 

doo = 1 — COS Ol — COS" 02 — COS 03 — COS 04 + (cos Ol COS 03 — COS 02 COS 04) 

The same holds for the configuration space of hyperbolic quadrilaterals, 
with cos and sin replaced by cosh and sinh. 

Proof. Consider the Gram matrix of the vectors from the center of the sphere 
(or center of the hyperboloid, in the hyperbolic case) to the vertices of the 
quadrilateral. Its determinant must be equal to zero. □ 

In the hyperbolic case, a similar equation is derived by considering the 
Minkowski space model of the hyperbolic plane. 

Theorem 7 (Ivory’s theorem on the sphere and in the hyperbolic plane). 

Let Ai and A 2 be two points on the sphere, respectively in the hyperbolic 
plane. Define an ellipse with the foci A\ and A 2 as the set of points with a 
constant sum of geodesic distances to A] and A 2 , and a hyperbola as the set 
of points with a constant absolute value of the difference of these distances. 
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Then in every curvilinear quadrilateral cut out by a confocal family of 
ellipses and hyperbolas the diagonals have equal geodesic lengths. 

Note that a spherical hyperbola consists of two spherical ellipses (replace 
one of the foci by its antipode). 


Proof. Coefficients of the equation ( 7 . 2 ) don’t change when a* are replaced 


by a t . For the coefficients dn, d\o, doi this is immediate from Lemma 7.1 
while doo must first be shown to be equal to 


— 1 + (sin a\ sin 03 + sin a 2 sin 04) 2 

— 2 sin a\ sin gq sin 03 sin aq — 2 cos oq cos a 2 cos 03 cos 04 

(with hyperbolic trigonometric functions in the hyperbolic case). Now the 
Ivory’s theorem follows by the same argument as in the euclidean case, see 
Section 16.11 □ 


It is not clear to us whether Theorem [ 7 ] is related to a version of Ivory’s 
theorem for hyperbolic spaces presented in [ 28 ]. 

Finally, note that while in the euclidean case the scaling of the sides 

(aq, 02, 03, 04) 1 y (Aai, A02, A03, A04) 

doesn’t change the configuration space, in the spherical and hyperbolic cases 
it does. 


8 . Open questions 

Problem 1 . The area S is a meromorphic function on the space of quadri¬ 
laterals with fixed side lengths. The inscribed quadrilateral 1 (a) is a branch 
point of A, and we have S(I(a )) = ^01020304 (Brahmagupta’s formula). 
What are the other branch points and other critical values of A 7 Is there a 
nice theory in the spherical and hyperbolic case? 

Problem 2 . Derive periodicity conditions for foldings of spherical and hy¬ 
perbolic quadrilaterals. 

Problem 3 . Study flexible bipyramids (or suspensions) with the help of 
parametrizations of the configuration spaces of spherical quadrilaterals. 

A bipyramid over a quadrangle is an octahedron, thus this is a gener¬ 
alization of Bricard’s problem of flexible octahedra. Flexible bipyramids 
were studied in ®E 3 . It is not known whether there exist non-trivial flexi¬ 
ble bipyramids over the pentagon; in the articles cited examples of flexible 
bipyramids over the hexagon were found. By cutting off the two apices of 
a bipyramid one obtains a polyhedron with vertices of degree four. This 
leads to a system of biquadratic equations, see Section [ 673 } Besides this sys¬ 
tem, there is also a holonomy condition (since the polyhedron is not simply 
connected). 
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Problem 4 . The dihedral group action on a = (01,02,03,04) does not 
change the configuration space Q^ r (a). What happens when the cyclic order 
of sides is changed in one of the following two ways? 

a' = (04,03,02,04), a" = (01,02,04,03) 



Figure 16 . Linking the configuration spaces Q^ r (a), 
Qor( a ')i and Qor( a ")■ 


The curves Q^ r (a), Q^ r (a'), and Qg r (a") are isomorphic, as follows from 
Theorem [ 2 ] or as is indicated on Figure 16 What does change, are the 


parametrizations. As formulas in Proposition | 3 . 12 ] suggest, the shifts start 
mixing up with the amplitudes. The configuration space of the linkage on 


Figure 16 seems to be a four-fold covering of each of the spaces Qff r (a), 


Qor(a'), and Q„ r (o" 
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